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Abstract 

We study the production of photons accompanied by jets in large-Q 2 deep inelastic 
scattering. Numerical results for the cross section differential with respect to the 
fraction of momentum z 7 carried by a photon inside a jet at large z 7 , up to 0(aa s ) in 
perturbative QCD, are presented. The sensitivity to the fragmentation contribution 
allows one to study the quark-to-photon fragmentation function. Our results can 
be confronted with future experimental data from HERA. 
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1 Introduction 



The production of final state photons at large transverse momenta in high-energy pro- 
cesses is an important observable for testing QCD. Data on large-p-r photon production 
in hadronic collisions have been used in the past to obtain information on the gluon dis- 
tribution in the photon. In addition, a good understanding of direct photon production 
within the framework of the standard model will be essential for new physics searches at 
future colliders. 

Photons from high-energy collisions are produced essentially by two mechanisms: the 
direct production of a photon off a primary quark or antiquark, or the fragmentation of 
a hadronic jet into a single photon carrying a large fraction of the jet energy. The direct 
production gives rise to perturbatively calculable short-distance contributions whereas 
the fragmentation contribution is due primarily to a long distance process which can 
not be calculated completely inside perturbative QCD. The latter is described by the 
process-independent quark-, antiquark-, or gluon-to-photon fragmentation functions (FF) 
HI which must be fixed by experimental data as long as they can not be calculated by non- 
perturbative methods. Their evolution with the factorization scale /ip and their so-called 
point-like contribution (up to a normalization) can, however, be calculated perturbatively. 

Directly produced photons are usually well separated from hadron jets, while photons 
originating from the fragmentation process are primarily found inside hadronic jets. By 
imposing an isolation criterion on the photon, one can, in principle, suppress (but in 
general, not eliminate) the fragmentation contribution to the final state photon cross 
sectionP]. 

So far, only a limited number of measurements of single photon production exists through 
which direct information on the photon FF can be obtained. A possible way is the 
measurement of the inclusive photon cross section in e + e~ annihilation. Recently the 
OPAL collaboration at LEP has measured this cross section for final state photons in the 
range 0.2 < x 1 < 1.0, where x 7 = 2E 1 /Mz is the photon energy fraction in terms of 
the beam energy ||. The results are in reasonable agreement with predictions obtained 
within perturbative QCD using the model-dependent parametrizations of the photon FF 
of DO U, GRV [| or BFG [|, when choosing the factorization scale fi F = M z . The DO 
model for the FF is based on an asymptotic solution of the evolution equation whereas the 
other two models H Q contain in addition a non-perturbative input inspired by the vector 
meson dominance model. Unfortunately, the experimental precision was not sufficiently 
high to discriminate between these theoretical predictions. Note also that in ]7|, |8| it 
was shown that predictions obtained within perturbative QCD and using a fixed-order 
expanded expression for the photon FF also agreed well with the OPAL data. 

An alternative way to determine the FF from e + e~ data is to measure the production of 



photons accompanied with a definite number of hadronic jets [10|. In this approach, the 



photon is treated like any other hadron and is clustered simultaneously with the other 



L See (|] for a proposal to eliminate the dependence on the photon FF. 



2 



hadrons into jets (so-called democratic clustering procedure) |TT[. Then, one of the jets in 
the final state is considered as the photon jet if the fraction of the electromagnetic energy 
inside the jet is sufficiently large, i.e. 



E 

z = — 7 > z cut (1) 

-£^7 + n had 

with z cu t fixed by the experimental conditions. On this basis the ALEPH collaboration 



at LEP analysed events produced at the Z-resonance |L2| which contained one hadron jet 
and a jet with a photon carrying at least 70% of the jet energy. This 7 + 1-jet rate was 
used to determine the quark-to-photon FF, which was calculated in leading order (LO), 
i.e. up to 0(a), in fllTf . The calculation was extended to next-to-leading order (NLO), i.e. 
up to 0(aa s ), in M, |J and a NLO photon-to-quark FF was obtained [0] by comparing 
with the ALEPH data. Due to the unfavorable signal-to-background ratio below z = 0.7 
for this particular observable, the 7 + 1-jet rate, the photon FF could be determined only 



in the range 0.7 < z < 1.0. In [§, [13|] also the predictions of GRV and BFG for 



the photon FF were compared to the ALEPH 7 + 1-jet data. It turned out that the BFG 
parametrization is in agreement with the measurements while the rate predicted with the 
help of the FF of GRV lies too high. 

The photon FF is process independent as any other FF and can be used to predict the cross 
section for the production of single photons in other processes. This has been done recently 
for isolated photon production in deep-inelastic ep scattering (DIS) jl4}| . In that work 
the emphasis was put on the prediction of cross sections in the HERA kinematic range 
without the necessity to introduce photon-parton separation cuts as used in earlier work 



[15]. This was achieved by including the divergent photon-quark (-antiquark) singular 
contributions and the bare quark-to-photon fragmentation contributions, leaving in the 
final result only a finite factorization scale dependent quark-to-photon FF contribution. 
In our previous work we adopted as the photon FF the lowest-order fit obtained from the 
photon+l-jet data of the ALEPH collaboration [|12| . 



In the present work we study the influence of the fragmentation contribution on the cross 
section for the production of a photon plus jets in DIS and give predictions for da/dz.y as 
a function of 27 in the interval z cut = 0.5 < z 1 < 1.0 (for the definition of 27 see below). 
For this purpose we calculate da/dz^ with the same cuts on the DIS variables as in [|14[] 
for essentially three different sets of photon FF's, the ALEPH sets in LO and NLO |12], |7]], 
the leading-logarithmic (LL) and beyond-leading logarithmic (BLL) parametrizations of 
GRV || and the BLL parametrization of BFG . The first set of photon FF is obtained 
within a fixed-order framework at a given order in a s . The two other sets were obtained 
after the leading and/or next-to- leading logarithms of the factorization scale fip were 
resummed. The essence of these two approaches has been described extensively in ]§[] and 
will be outlined below. 
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2 Leading and Next-to-Leading Order Cross Section 



In leading ordeiQ (i.e., at 0(a)), the cross section for the production of photon plus jets in 
DIS receives contributions from the quark (antiquark) subprocess 7*g — > jq (7*5 — > jq). 
Together with the remnant jet from the proton this process gives rise to a 7 + (1 + 1)- 
jet final state. In the virtual photon-proton center-of-mass system the hard 7 recoils 
against the hard jet back-to-back. Cuts on the usual DIS variables x, y and Q 2 are 
applied to remove 7 production by incoming 7*'s of small virtuality. To produce 7's 
of sufficiently high energy, an explicit cut on the total j*p center-of-mass energy, W, is 
introduced. Both leptons and quarks emit photons. The leptonic radiation is suppressed 
by a cut on the photon emission angle with respect to the incoming electron in the same 



way as in the earlier work [FJ, [TjJ. In LO each parton is identified with a jet and the 
photon is automatically isolated from the quark jet by demanding a non-zero transverse 
momentum of the photon or jet in the 7*p center-of-mass frame. At this order there is no 
fragmentation contribution. 

At next-to-leading order (0(aa s )) we have subprocesses with an additional gluon, either 
in the final or in the initial state, i.e. the subprocesses 7*g — > ^qg (and similarly with 
q replaced by q) and 7*^ — > •yqq, respectively. In addition, virtual corrections (one-loop 
diagrams at 0(aa s )) to the LO processes have to be included. The processes 7*g — > jqg 
and 7*g — > ^qq contribute both to the 7 + (1 + l)-jets cross section and to the cross 
section for 7 + (2 + l)-jets. In the latter case each parton in the final state (including the 
photon) builds a jet on its own, whereas for 7 + (1 + l)-jets a pair of final state partons 
is combined into one jet. The recombination of two partons will be performed with the 
help of the cone algorithm. The exact prescriptions will be given later when we present 
our results. In the NLO calculation of the 7 + (1 + l)-jet cross section one encounters 
the well known infrared and collinear singularities. For the processes 7*g — > 'yqg and 
1*9 ~~ > 1QQ they appear in those phase space regions where two partons are degenerate 
to one parton, i.e. when one of the partons becomes soft or two partons become collinear 
to each other. These singularities cancel in the case of soft gluons or for collinear quark- 
gluon pairs against the singularities from the virtual corrections to the LO process or 
have to be factorized and absorbed into the renormalized parton distribution functions 
(PDF's) of the proton. To accomplish this cancellation, the singularities are isolated in 
an analytic calculation with the help of dimensional regularization and the phase space 



slicing method. The technical details of this procedure have been described in p4| , [L5 
and need not to be repeated here. 

The matrix elements \M\ 2 for the processes 7*g — > 'yqg and j*g — * 'jqq have also photonic 
infrared and collinear singularities. The infrared singularity is outside the kinematical 
region we are interested in since the photon is required to be observed in the detector. In 
the numerical calculation of the cross section we impose this condition by demanding a 
minimal transverse momentum of the final state photon. This cut removes also collinear 

2 Here and in the following, we do not count the extra factor a from the e — 7* vertex. 
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singularities due to initial state photon radiation. Final state collinear singularities are 
present and are treated with the phase space slicing method similar to the case of the 
quark-gluon collinear contribution. For this purpose the phase space slicing parameter 
Uq is introduced, which is chosen very small, so that the matrix element \M\ 2 can be 
approximated by its singular part. The slicing of the phase space is done with the help 
of the scaled squared invariant masses, for example, with y iq = (jp q + p<y) 2 /W 2 for the 
subprocess 7*g — > jqg, where p 7 and p q are the four-momenta of the outgoing photon 
and quark, respectively. In the gluon initiated subprocess one has two singular regions, 
which are controlled by the variables y iq and y iqi respectively. In the region y iq > yl 
the cross section is evaluated numerically. The details for the calculation of the various 
contributions are described in [14| . 



The contribution to \M\ 2 in the region y iq < is collinearly divergent and is regulated 
by dimensional regularization. The divergent part is absorbed into the bare photon FF 
to yield the renormalized photon FF denoted by D q _+ T The additional fragmentation 
contribution to \M\ 2 for the subprocess 7*g — > ^qg has the following form 

\M\ 2 rq ^ qg = \M\ 2 rq ^ qg ®D q ^{z). (2) 

There is a similar expression for the subprocess 7*g — > ^qq. It is obvious from (Q) that 
the fragmentation contribution is 0(aa s ), the photon FF D q ^(z) given by 

D q ^(z) = D q ^(z, &) + ^ (p£\z) In Z{1 ~ Z \ ylW2 + z] (3) 

being of 0(a). D q ^(z, fMp) in (Q) stands for the non-perturbative FF of the transition 
q — > 7 at the factorization scale fip, i.e. the scale at which the redefinition has been 
performed. This function will be given by one of the sets of photon FF mentioned above. 
The second term in (^), if substituted in (|2]), is the finite part resulting from adding the 
bare photon FF and the collinear photon-quark (-antiquark) contribution to the matrix 
element |M|L integrated in the region y iq < y^. 

The yl dependence in (||) is expected to cancel the dependence of the numerically eval- 
uated 7 + (1 + l)-jet cross section restricted to the region y iq > y^ as studied in |15 . 



This means that y$ is only a technical cut separating divergent from finite contributions. 
This collinear photon-quark (-antiquark) contribution to the matrix element |M| 2 „ g _^ 7(?5 
being calculated in the collinear approximation is valid only up to terms 0{yo). Con- 
sequently the cut ?/o must be chosen sufficiently small. The independence must be 
checked by varying t/J. Results of this test will be shown below. In @, P@>(z) is the LO 
quark-to-photon splitting function 

P£?{z) = 1 + (1 ~* )2 (4) 

H1 z 

and e q is the electric charge of quark q. The variable z denotes the fraction of the final 
photon energy in terms of the energy of the quark emitting the photon. If the photon is 
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emitted from a final state quark with four-momentum p q > = p qi + p 7 , then z is obtained 
from p 1 = zp q i 4 . It can also be related to the invariants y q3l and y q3q4 obtained from the 
four-momenta of the subprocess 7*^3 — > ^q^g with the result 

According to (0) the fragmentation contribution to the subprocess 7*^3 — > ^q^g is calcu- 
lated from the convolution of the FF with the 0(a s ) matrix element \M\ 2 of the process 
7*g 3 — > q±g, which is well known. It yields a term of 0(aa s ) to the cross section for 
7 + (1 + l)-jets. Similarly the fragmentation contribution to the subprocess 7*$ — > ^qq is 
calculated from the convolution of the FF for q — > 7 with the matrix element \M\ 2 of the 
process 7*^ — > gg which is also known. Equivalent formulas are used for the calculation 
of the fragmentation contribution to 7*^3 — > ^q^g and 7*g — > 7gg, where the g fragments 
into a photon. The fragmentation is always assumed to occur collinearly, i.e. with no 
additional transverse momentum in the fragmentation process. 



3 Quark-to-Photon Fragmentation Functions 

As long as we intend to calculate the 7 + (n + l)-jet cross section only up to 0(aa s ) we 
need the photon FF only in LO, i.e. at 0(a). This means that the gluon-to-photon FF, 
which starts at 0(aa s ), is not needed. In LO the non-perturbative quark-to-photon FF 
obeys the evolution equation 



dD q ^(z, /ip) _ ae| 
dln(i F 2tt 



PS}{z) (6) 



with the solution 



Z7T \ PO / 



(z,^ ). (7) 



D q ^ y (z, fi ) is an initial value which must be fitted to experimental data with a chosen 
initial scale [Iq. This has been done by the ALEPH collaboration From the fit to 
their e + e~ — > 7 + 1-jet data [0 they obtained 
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D^(z,^) = ^ (-Pff(z) ln(l - zf - 13.26 

with /i = 0.14 GeV. (|7|) together with (j8|) is one of the photon FF choices which will be 
used to predict the cross section da/dz^ for 7+ (n+l)-jet production in DIS ep scattering. 

We note that (|7|) is an exact solution of (||) at 0(a). Furthermore when we substitute (|7|) 
into (0) we see that together with the finite contribution (i.e., the second term in @) the 
cross section becomes independent of the factorization scale fiF. This means that for the 
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cancellation of the fj, F dependence only the LO FF is needed. Nonetheless, in order to see 
the influence of the NLO corrections to D q ^(z, fi F ) we shall evaluate the 7 + (n + l)-jet 
cross section also with the inclusion of the NLO photon FF. 

Similarly to ([?[) the NLO FF D q ^(z, fi F ) is obtained as the solution of ([]), but now with 
0(aa s ) terms added on the right-hand side of (||). The result at scale fi F is 



D 



a 



ae 



2 r 



2tt 




>>Mo) 



(9) 



a 



where P^'(z) is the next-to-leading order quark-to-photon splitting function [16] and 



Pj?'(z) is the well-known LO qq splitting function. D q ^(z, /io) is the initial value of the 
NLO FF, which contains all unknown long-distance contributions. The result in (J|) is an 
exact solution of the evolution equation up to 0(aa s ). Based on the above arguments, 
also the NLO photon FF has recently been determined using the ALEPH 7 + 1-jet data 
||12|| . A three parameter fit with a s = 0.124 (this value for a s , just a scale-independent 
parameter here, was chosen so that the observed total e + e~ annihilation cross section into 
hadrons is reproduced in the 0(a s ) calculation) yields M 



2tt 



-P$(z) ln(l - z) 2 + 20.8(1 -z)- 11.07 



(10) 



with fiQ = 0.64 GeV. Inside the experimental errors this fit for the photon FF at /i 
describes the ALEPH data as good as the LO fit (|]) [7| [| . A fit with a larger value of a s 
is also reported in ||. 

It should be noted that the LO and NLO FF's of the photon as given in ([?[) together 
with ([D and in @ together with (0) do not take into account resummation of powers of 
ln(/i|i//io) as usually implemented, e.g. via Altarelli-Parisi evolution [r?|. Such resumma- 



tions are only unambiguous if the resummed logarithm is the only large logarithm in the 
kinematical region under consideration. If logarithms of different arguments can become 
simultaneously large, the resummation of one of these logarithms at a given order implies 
that all other potentially large logarithms are shifted into a higher order of the pertur- 
bative expansion, i.e. are neglected. In the evaluation of the 7 + 1-jet rate at 0(a) 
and 0(aa s ) [T]J at LEP for 0.7 < z < 1, one encounters at least two different potentially 
large logarithms, In fi F and ln(l — z). In the high z region, the region where the photon 
is isolated or almost isolated it is by far not clear that ln/x^ is the largest logarithm. 
Choosing not to resum the logarithms of Infip is therefore equally justified for the case 
of large z — > 1. 

In the conventional approach, the parton-to-photon FF's Di^(z, fi F ) satisfy a system 



of inhomogeneous evolution equations |T7| . Usually these equations are diagonalized in 
terms of the singlet and non-singlet quark FF's as well as the gluon FF. For the case that 
the gluon-to-photon fragmentation is neglected, which we shall do in our application to 
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DIS 7+jets production, these equations can be simplified J8|. Then the flavor singlet and 
non-singlet quark-to-photon FF's become equal to a unique function -D 9 _> 7 . This function 
satisfies the all-orders evolution equation 



which has a similar form as the next-to-leading order evolution in the fixed-order approach, 
but the coupling a s is not fixed and now is taken as a function of the factorization scale 
/jlf- The full solution D q ^ 1 (z, fj,p) of the all-orders inhomogeneous evolution equation 
is a sum of two contributions. The first term is the point-like (or perturbative) part 
D? 1 (z, fip), which in the leading logarithmic (LL) approximation is a solution of the 
inhomogeneous equation (|TT|) . The second term is the hadronic (or non-perturbative) 
part (z, fip), which is the solution of the corresponding homogeneous equation. This 
term must be fitted to experimental data to obtain the input at a starting scale. At LL 
only terms of the form a" ln n+1 jx F are kept while beyond leading logarithms (BLL) both 
the leading and the sub- leading (a"ln n /ip) logarithms of fip are resummed to all orders 
in the strong coupling constant a s . Thus the LO solution ([/]) is the n = term of the LL 
approximation, whereas the NLO solution (|9]) is the BLL approximation up to the order 
n = 1. By dropping the a s terms in (H) one can also infer the LL approximation up 
to the order n — 1. 

In this conventional approach, the quark-to-photon fragmentation function is regarded as 
being of 0(a/a s ). In fact, the n = term of the LL approximation is proportional to 
\n(fi 2 F /fil). At large scales [xp 3> /ig and assuming that fip can be identified with the 
scale that determines a s , i.e. a s ~ l/ln/z^, the LL FF becomes of the order of 0(a/a s ). 
This motivates the usual statement that the total fragmentation contribution is like 0(a), 
i.e. of the same order as the LO direct photon contribution. Of course, this statement 
depends on the considered kinematic region. With this argument however, one usually 
justifies the inclusion of terms in the fragmentation contribution of 0{aa 2 s ) which become 
of 0(aa s ) if the FF is considered of 0(a/a s ). 

The most recent parametrizations of photon FF's in the conventional approach as de- 
scribed in the previous paragraph are those of GRV || and BFG ||. The FF's of GRV 
H exist in the LL as well as in the BLL form, which in the following we shall denote as 
the LO and NLO FF's of GRV. However, the FF of BFG § are available only in the BLL 



approximation. They both have been compared to the ALEPH 7 + 1-jet cross section fl2fl , 
which is sensitive to the large z region (0.7 < z < 1.0). The BFG prediction goes through 
the experimental points whereas the GRV NLO parametrization lies systematically higher 
and agrees less well with the ALEPH data || [TB|. This difference can be attributed to 
both the choice of the input scale /io and of the non-perturbative input at this scale. The 
BFG input is smaller and according to the ALEPH data this choice is preferred. 

Thus we have at least five different versions of photon FF's at our disposal which have 
been compared to the ALEPH data: two LO parametrizations, the one using the ALEPH 
data directly to determine the initial distribution given in (f7|) and (||) and the GRV 
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parametrization in LO; and three NLO choices, the one written in (^j and (|T0| ) directly 
fitted to the ALEPH data, and the NLO parametrizations of GRV and BFG. We shall 
use these five parametrizations to predict the 7 + (n + l)-jet cross section in deep-inelastic 
ep scattering. 



4 Numerical Results 

The results for these cross sections are obtained for energies and kinematical cuts appro- 
priate for the HERA experiments. The energies of the incoming electron and proton are 
E e = 27.5 GeV and E p = 820 GeV, respectively. The cuts on the DIS variables are as in 



our previous work fl4| 



Q 2 >10GeV 2 , H/>10GeV, 
1(T 4 < x < 0.5, 0.05 < y < 0.95 . (12) 



To reduce the background from lepton radiation [15] we require 



90° < B 1 < 173°, 9 ie > 10° (13) 

where 9 1 is the emission angle of the photon measured with respect to the momentum of 
the incoming electron in the HERA laboratory frame. The cut on 9 iei the angle between 
the photon and the outgoing electron momentum, suppresses radiation from the final-state 
electron. The 7 and the hadron jet J are required to have minimal transverse momenta 

VT n > 5 GeV, p T:J > 6 GeV . (14) 

Different values of minimal pr's for the photon and the jet have to be chosen in order to 
avoid the otherwise present infrared sensitivity of the NLO predictions ffM. The PDF's 



of the proton are taken from MRST ||18|| . a s is calculated from the two- loop formula with 
the same A value (Aj^g(nf = 4) = 300 MeV) as used in the MRST parametrization. The 
scale in a s and the factorization scales fip for the proton PDF and the photon FF are 
equal and fixed to y/Q^. 

We are interested in the differential cross section da/dz.y at NLO (up to 0(aa s )) as a 
function of z 7 , where 

z, = ^^ . (15) 

VT n + PT,had 

This definition of z 1 agrees with the one in (||) for the fragmentation contribution. In ([15]) , 
Pr,had is the transverse momentum of the parton producing hadrons, which is recombined 
with the photon into the photon jet. Here the photon is treated like any other parton 
during recombination, so that one of the recombined jets may be the photon jet. For the 
recombination of the two partons into a hadron jet or a parton and the photon into the 



photon jet, we use the cone algorithm of the Snowmass convention [19|. The recombination 
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Figure 1: yj-dependence of separate contributions to the 7 + (n + l)-jets cross section. 

is applied in the 7 + 2 parton sample and yields a contribution to the 7 + (1 + l)-jet class. 
The 7 + (2 + l)-jet class consists of the unrecombined contributions of the 7 + 2 parton 
sample. The recombination is performed in the 7*p center-of-mass frame. Two partons i 
and j are combined into a jet J if they obey the cone constrains R it j < R and R^j < R, 
where 

Ri,j= v/(^-^) 2 + (0,-0j) 2 . (16) 

r/j and j are the rapidity and the azimuthal angle of the recombined jet. If for example, 
parton % is the 7, then J is the photon jet. The jet variables rjj and <pj are obtained from 
the averages of the corresponding variables of the recombined partons % and j multiplied 
with their respective px values and pt,j = Vt,% +Ptj- We choose R — 1. The azimuthal 
angle is defined with respect to the scattering plane defined by the momenta of the ingoing 
and outgoing electron. It is known that the cone algorithm is ambiguous for final states 
with more than three particles or partons. Since we have maximally three partons in 
the final state this is not relevant in our case. Furthermore it will be no problem to 
repeat the calculation for any other cluster algorithm that might be used in the analysis 
of forthcoming experimental data. It is obvious that for the fragmentation contribution 
the cone constraint is always satisfied since the hadronic remnant in the fragmentation is 
collinear with the photon. 

Before we present our results for da/dzy for the various choices of the photon FF we 
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da(-y + (n + 1)- jets) / dzy [pb] 




Figure 2: Comparison of the 7 + (n + l)-jet cross section for the LO and NLO parametriza- 
tions of the FF obtained from fits to the ALEPH data. 



show that the cross section is independent of the slicing cut y^. For this purpose we have 
calculated the sum of the cross sections for the production of 7 + (n + l)-jets (n = 1, 2) 
as a function of yj in the range 2 • 10~ 7 < y$ < 10~ 3 . The result is shown in Fig. |I|. Here 
we have plotted the cross section for the direct contribution with y 7 > y$ which increases 
with decreasing y$ and the direct contribution for y 1 < together with the second term 
in contribution (|3]) inserted in (0) which we denote by 'frag(sing)' in the figure. This latter 
contribution is negative and decreases with decreasing y$ in such a way that the sum of 
the increasing and decreasing contributions is constant inside the numerical accuracy for 
< 5 ■ 10~ 4 . The sum shown in Fig. H contains also the fragmentation contribution @ 
with the photon FF D q - n (z J fi%) as given in ([?[) and (|J) (denoted by 'frag(fin)' in the 
figure). Of course, this contribution is independent of y$. It amounts to approximately 
14% of the total 7 + (n + l)-jet cross section. The cross sections shown in Fig. [I] are 
calculated with the kinematical cuts specified above and the additional cut z 1 > 0.5. 
From Fig. [I] we conclude that the slicing cut for the photon should be chosen smaller than 
10~ 4 to obtain a reliable cross section. For the following results we fixed it to = 10~ 5 
as in our previous work [13 . 

The results for da jdz 1 as a function of z 1 are shown in Figs. 0, ^| and (|. The total cross 
section is dominated by the last bin 0.975 < z 7 < 1 which is clipped off above 10 pb 
in these figures. In this bin we have the direct contribution and those fragmentation 
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contributions with no or little hadronic remnant. The fraction of the total cross section 
(for 0.5 < z 7 < 1) contained in the last bin amounts to 70 % to 80 % and does not depend 
on the parametrization of the FF. Therefore it is clear that the dependence of the photon 
plus jet cross section on z 7 below z 7 = 1 is an appropriate observable which contains 
information on the photon FF. Here the fragmentation contribution dominates and in 
addition we have those direct contributions from the (7 + 2)-parton sample where the 
photon is recombined with one of the partons. In Fig. |2| the cross section da /dz 1 in the 
region 0.5 < z 1 < 1.0 is plotted for the photon FF's from (0) and (||) (denoted ALEPH) 
and from (f|) and fllOD (denoted GGdR(a s = 0.124) in the figure). For both cases the cross 
section first decreases with increasing z 7 until it reaches a minimum near z 7 = 0.9, from 
which it increases strongly towards z 1 = 1. This dependence is similar to what has been 
found for the photon plus one-jet cross section in e + e~ annihilation at the Z-resonance 
|0, p. Note however, that below z 1 = 0.5, i.e. outside the range shown in these figures, the 
cross section decreases again towards smaller values of z 1 . Unlike at LEP, the photon-jet 
is here required to have a minimal transverse momentum as given by ([15|). The decrease 
in the cross section for values of z 1 below 0.5 is a kinematical effect resulting from the 
imposition of this cut. 

For the two choices of FF's the predicted cross section differs by approximately 20%. The 
cross section for the LO FF is smaller than for the NLO FF. We note that the a s value 
quoted for this FF is not the value for which the cross section was calculated; rather it 
is the value which was chosen to fit the photon FF to the ALEPH data. Both photon 
FF's fit the ALEPH data in LO and NLO approximation. In our case only the LO 
FF is appropriate. The result for the NLO FF gives just an indication of higher-order 
corrections up to 0(otot 2 s ). However, it is not a full prediction at this order since the 
NNLO calculation for the direct contribution has not been done, the NLO corrections to 
\M\ 2 on the right-hand-side of (|2|) are not included either. 

In Fig. [3] the predictions with GRV |5|] and BFG || FF's are shown. Both are FF's in 
NLO and with all leading and subleading logarithms of \\i([i 2 F /mfy resummed at the next- 
to-leading logarithmic accuracy. The behavior of the cross section as a function of z 1 is 
qualitatively similar to the result in Fig. |2]. The cross section for GRV is larger than the 
one with the BFG FF. This is to be expected since GRV gives a larger e + e~ cross section 
for (7 + l)-jet than BFG also for the ALEPH kinematical conditions. 

Comparing with the fixed-order result in Fig. |2| the BFG cross section lies above the fixed- 
order prediction obtained using the LO FF over the whole z 7 region. The discrepancy is 
largest in the region 0.65 < 2 7 < 0.95 where the cross section is smaller. At the minimum 
the two cross sections differ by more than a factor two. 

In Fig. [| the results for the GRV photon FF's are displayed, but now for both the LO 
and the NLO parametrization. On the average the cross section for the LO fragmentation 
function is again smaller by 20% compared to the NLO prediction. In addition, we have 
plotted the cross sections for the point-like approximations of the FF's. This approxi- 
mation is rather good, i.e. the influence of the hadronic part of the FF is small in the 
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Figure 3: Comparison of the 7 + (n + l)-jet cross section for the NLO parametrizations 
of GRV and BFG. 



considered z 7 range. We note that the difference between the full FF and the point-like 
part increases with decreasing z 1 . 

If we compare our results for the various photon FF's in Figs. 0, ^| and [| we observe that 
the predictions agree approximately within 20% in the small and large z 1 regions, i.e. for 
z 7 < 0.65 and z 7 > 0.90. However, near the minimum of the cross section, i.e. in the region 
0.65 < £ 7 < 0.95, the results differ by up to a factor two. The largest differences occur 
between the predictions obtained with the ALEPH photon fragmentation function on the 
one side and the GRV and BFG parametrizations on the other side. This difference comes 
mainly from the fact that different evolution approaches are used. Whereas for GRV and 
BFG the FF at fi F = Q 2 is obtained from the conventional evolution after the leading 
and/or subleading logarithms of \xp were resummed, the ALEPH FF's are evolved only to 
the respective finite order in a s as given in (^) and (|9|). Therefore, if we calculated da/dz^ 
for ep — > e7 + (n + l)-jets at the large scale \lf — Mz the cross sections obtained for 
BFG and ALEPH would come out quite similar over the whole 27 range inside the 20% 
margin. Only when we go to the scale \j? F = Q 2 , which on the average is much smaller, 
we observe that the cross section obtained using the BFG photon fragmentation function 
is much larger than the ALEPH cross section in the region 0.65 < z 1 < 0.95. 

Provided the resummed solution of the all-orders evolution equation can be accurately 
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Figure 4: Comparison of the 7 + (n+l)-jet cross section for LO and NLO parametrizations 
of GRV. 



determined || over the whole z 1 range under consideration, i.e. for 0.5 < z y < 1, the 
approach using this solution represents the theoretically preferred approach as it is the 
most complete. The approach using an expanded and therefore approximated photon FF 
has however important advantages. As already mentioned, its use leads to factorization 
scale independent results for the cross section evaluated at a given fixed order in a s . 
Moreover it enables an analytic determination of the photon FF. 

Within the conventional approaches, the BFG prediction should be preferred over the 
result with the GRV FF. The difference between GRV and BFG in Fig. is related to 
the choice of a different input at the starting scale /j,q. As already mentioned, the BFG 
parametrization is in better agreement with the ALEPH data || than the GRV NLO 
parametrization and therefore is the more realistic choice of FF for the photon. This 
could be further tested as soon as data from HERA on photon plus jet production in DIS 
become available. 
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5 Conclusions 

We have presented a NLO calculation for the production of photons accompanied by jets in 
deep inelastic electron proton scattering, taking into account the contribution from quark- 
to-photon fragmentation. We have calculated the cross section as a function of z 1 for those 
parametrizations of the photon fragmentation functions, which had been compared to 
photon plus jet data from the ALEPH collaboration. As we observed significant differences 
between the predictions, we expect that the measurement of photon plus jet production 
in DIS at HERA will contribute to testing these photon fragmentation functions. 
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